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We have studied the time evolution of magnetization of a high spin magnetic cluster in the ground 
state, in the presence of a sinusoidal axial magnetic field and a static transverse field by explic- 
itly solving time dependent schrodinger equation. We observe oscillations of magnetization in the 
plateau region which has all the characteristics reminiscent of the oscillation of probability of occu- 
pation of a state in a two level lattice in the presence of an oscillating electric field. The high spin 
of the ground state leads to a large but finite two-level pseudo-lattice. The oscillations in magneti- 
zation occur at amplitude of magnetic fields achievable in a laboratory and they also persist for a 
wide range of spin-dipolar interactions. Thus, our study provides a magnetic analog of the optical 
two-level lattice. 



PACS numbers: 75.45.+j, 75.20.-g, 75.50.Xx, 75.40.Gb, 42.50.Gy 



In recent years, synthesis of high nuclearity transition 
metal complexes in high-spin ground state has spurred 
interest in magnetism on a nanoscale fjj]. The synthesis 
of Mni2 and Feg clusters in S=10 ground state has led 
to extensive study of quantum resonant tunneling and 
quantum interference phenomena The quantum res- 
onant tunneling manifests as plateaus in the magnetiza- 
tion vs magnetic field plots, with the width and location 
of plateaus being determined by the ramping speed of the 
magnetic field as well as the initial state ||. The quan- 
tum interference phenomena observed in the Fes clus- 
ter is because the paths connecting the Ms = +10 and 
Ms = — 10 could interfere in the presence of a magnetic 
field, leading to an oscillation in the tunneling probabil- 
ities j§. 

In these systems, the magnetic interaction between 
transition metal ions in neighbouring molecules is weak 
due to the bulky nature of the molecular complex. Thus, 
we can treat these molecules as isolated magnetic clus- 
ters. These clusters are characterized by exchange cou- 
plings between magnetic ions that are frustrated, leading 
to a high spin ground state separated from other low- 
lying high spin states by rather small energy gaps. These 
magnetic clusters also exhibit signatures of spin dipolar 
and other higher order spin-spin interactions allowed by 
the symmetry of the cluster. 

The exchange Hamiltonian of the Mni2 system with 
which we shall be concerned in this letter can be written 
as 

H cxch = Jij§i ■ Sj (1) 

<ij> 

where the exchange interactions are given in fig. 1 . The 
model exact S=10 ground state as well as the low-lying 
excited states of this Hamiltonian have been obtained 
from this model by using all the spatial symmetries and 



partial spin symmetry adaptation. The ground state and 
the first few low-lying states live in different symmetry 
subspaces and do not mix unless the molecular symmetry 
is broken H. Therefore, for computing properties such 
as quantum resonant tunneling at low-temperatures, it is 
sufficient to work with the different Ms manifolds within 
the S = 10 ground state. The spin-spin interactions al- 
lowed by the symmetry of the molecule leads to a spin 
Hamiltonian in the presence of an applied magnetic field 
given by 

Hdip = D S^ total + c (S^ total + Sy total ) 

~t~ Scrown H(t ) • S cr0 wn ~t~ <?corc H(t ) ■ S corc (2) 

where, for chemical reasons, the g factors for the Mn 3+ 
and Mn 4+ ions in the crown and core of the molecule (fig 
1) are slightly different A time-dependent solution of 
the above Hamiltonian, starting from an initial state with 
S = 10 and Ms=-10 shows quantum resonant tunneling 
when the magnetic field is ramped. || 

The energy level structure of the eigenstates of eqn. 
2 are very similar to those encountered in the dynamic 
studies of two-level lattices. There has been considerable 
study of the two-level lattices in which Stark- Wannier 
effect and Bloch oscillations are observed In the 

presence of oscillatory electric fields, population trapping 
and dynamic localization [|l0| have also been observed. 
The aim of the present study is to see if these effects 
can also be observed in the Mni2 magnetic cluster in the 
presence of an oscillatory magnetic field. 

In the present work we have studied the magnetiza- 
tion oscillations under the influence of a strong axial ac 
magnetic field and a weak stationary transverse field, us- 
ing explicit time evolution of an initial state. We start 
with an initial state which is the lowest eigenstate of the 
full exchange Hamiltonian (eqn. 1) with S = 10 and 
A/5 = —10. In the absence of the quartic spin terms and 
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external transverse field, this is also the lowest energy 
state of the low-energy Hamiltonian (eqn. 2). We evolve 
this state with time explicitly, using the time-dependent 
schrodinger equation, 

ih^=H dip .V>(t). (3) 

At each time step, we calculate the magnetization of the 
state. The value of D estimated from experiments is 
0.56K, while c is 3 x 10~ 5 K ||. The g cr0 wn and g COTC are 
respectively 1.96 and 2.00. In our studies, we have var- 
ied the amplitude of the time dependent magnetic field 
between 10 to 100 D. These fields range from 3.75 T to 
37.5 T , for experimentally estimated D value in Mni2 , 
which are indeed attainable The transverse field is 
static and fixed at a value equal to one hundredth the 
amplitude of the time varying axial magnetic field. 

The initial state in the time evolution is chosen to be 
the state with S — 10 and Ms = -10, which is the ground 
state, in the absence of the weak off-diagonal terms of 
Hdip. The angular frequency, u>, of the axial field is var- 
ied between 10 _1 and 10~ 3 radian — D/ti. The time evo- 
lution is carried out successively in steps At = 0.1 and 
the evolution is carried out over several periods of the 
applied field. 

In fig. 2a, we show a plot of magnetization vs time, 
for the amplitude of the axial field, Ho=50D, and angu- 
lar frequency uj = 10~ 2 radian — D/h. We notice that 
the magnetization shows distinct plateaus and in each 
plateau, the magnetization oscillates a fixed number of 
times. The number of oscillations N osc . in each plateau 
is given by the ratio of half the energy gap between the 
Ms = — 10 and the Ms = —9 states at the field Ho and 
Huj, where w corresponds to the angular frequency of the 
applied axial field. This kind of oscillations with simi- 
lar dependence of N osc . has been theoretically observed 
in the two-level systems by Rotvig et al jl2j and Ragha- 
van et al Such oscillations have also been seen in 

the Bose-Einstein condensates in a double- well trap [|l4j . 
Rotvig et al observe this in the context of a two-band 
semiconductor superlattice in an external electric field, 
while Raghavan et al observe it in a single band model 
in the presence of an electric field. In these calculations, 
the probability in a given state shows the temporal oscil- 
lations that we see here for magnetization. Indeed, the 
different Ms states in the magnetic cluster can be viewed 
as forming a large but finite lattice. The transverse mag- 
netic field couples the states at successive lattice sites, 
much as the transfer terms in the models of Rotvig et al 
and Raghavan et al. The time varying axial magnetic 
field corresponds to the applied oscillatory electric field 
of the two level systems. It is also worth noting that 
for the parameters that Raghavan et al use, the oscilla- 
tions die down for larger lattices. It appears that the size 
of the pseudo-lattice provided by the magnetic cluster is 
not large enough for the oscillations to die down for the 
realistic model parameters we have chosen. 



The oscillations in magnetic field that we observe are 
quite robust. We have observed these oscillations for 
other initial states with integral Ms values corresponding 
to the eigenstates of the Hdi p ., with off-diagonal elements 
set to zero (fig. 2 b-d). The dependence of the magneti- 
zation on the amplitude of the magnetic field is shown in 
fig. 3. We note that at higher fields, we see more oscilla- 
tions in each plateau since the gap between successive Ms 
states widen with increasing amplitude. In fig. 4 b,d and 
f, we show the energies of the two low-lying eigenstates 
of Hdip., corresponding to the axial field at that instance. 
We note that the jump in magnetization coincides with 
two states with different Ms values becoming degenerate. 
In the time between these coincidences, the magnetiza- 
tion shows small amplitude oscillations. Assuming that 
the system wavefunction evolves as exp(-i(E2 — Ei)t/h), 
where Ei and E2 are the average energies of the two 
states in question, the number of oscillation in a time 
period t = 2tt/uj is given by (E2 — Ei)/hct;. We note that 
the energy gap at the maximum amplitude in all the cases 
are the same and hence the number of oscillations in a 
plateau is inversely proportional to the frequency of the 
axial field. In fig. 4 a,c and e we show the dependence 
of the oscillation pattern on the frequency of the axial 
field with a fixed amplitude. We note that at decreasing 
frequencies, the number of oscillations in each plateau 
increases and the plateau structure itself vanishes with 
the magnetization following the magnetic field for higher 
frequency of the axial field. 

We have also studied this system in the presence of 
an axial magnetic field with two different frequencies, 
given by H cos(wit)cos(w2t). This leads to a beat pat- 
tern involving the sum and difference of the two frequen- 
cies. In fig. 5, we show the time variation of magnetiza- 
tion in the presence of two different frequency magnetic 
fields. We note that one set of oscillations correspond to 
N osc =AE/2 h(uji +W2) while another set of oscillations 
have N osc .=Ai?/2 h(usi — u>2), where AE is the difference 
in energy between the Mg=-10 and Mg=-9 states when 
the axial field is equal to the amplitude. 

The parameters for which we have carried out the cal- 
culations corresponds to Mni2. However, systems such 
as Fes also have high spin ground state although the D 
value is different. In order to study the effect of change 
in D value on the oscillations, we have carried out these 
calculations for several Ho/D values, in fig. 6, we show 
the results of our calculations for Ho/D = 20 and Ho/D 
= 10. The oscillations vanish (Ho/D = 20) and we have 
aMret behaviour that does not have much structure for 
this larger Ho/D value. However for lower Ho/D value, 
we find that the oscillations persist, but with much re- 
duced amplitude. Thus, it is possible that in other high 
spin systems, these oscillations are found at different field 
amplitude. 

To conclude, the magnetic clusters in high spin ground 
state behave like a finite lattice of two level systems and 
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in the presence of an oscillating axial field and static 
transverse field lead to predictable oscillations in the 
plateau regions of magnetization vs time plots. The field 
at which these oscillations can be observed fall in the 
range of 4 to 40 T with a frequency ranging from 1 to 10 
MHz for the known high-spin magnetic clusters. While 
such behaviour is predicted in optical two level systems, 
their observation may be easier in high nuclearity high 
spin magnetic systems. 
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Figure Captions 



1 . A schematic diagram of the exchange interactions be- 
tween the Mn ions in the Mn^Ac molecule. The inter- 
actions J\ = 215K and J2 = J3 = 86K are antiferromag- 
netic, while J4 = — 64.5K is ferromagnetic. The Mn 4+ 
ions form the core while the Mn 3+ ions form the crown. 

i 

2. (a) Plot of evolution of magnetization with time (in 
units of h/ D x 10" 1 ) for H /D =50 and w=0.01 
radian — D/h, initial state is M s =-10. For (b), (c) and 
(d) Ho/D =100 and ut—0.01 radian — D/h, initial states 
are respectively M s =-9, -4 and 4. 

3. (a) Plot of magnetization vs time when Ho/D =100 
and w=0.01. For (b) H /D =200 and w=0.01. Initial 
state in both cases is M s =-10. w and time are in units 
of fig 2. 

4. Magnetization vs time (a, c and e) and energy of M s =- 
10 (solid line) and -9 (broken line) vs time (b, d and f) 
for (i) w=0.001 (a, b) (ii) w=0.01 (c, d) and (iii) w=0.1 
(e, f). Ho/D = 100 for all cases. u> and time are in units 
of fig 2. 

5. Plot of magnetization vs time when Ho/D =100, 
c<ji=0.01 and W2=0.005. to and time are in units of fig 
2. 

6. Evolution of magnetization with time with different 
H /D values, lu is kept fixed at 0.01. H /D=100, 20 and 
10 respectively for (a), (b) and (c). to and time are in 
units of fig 2. 
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